Andreev bound states and current-phase relations in three-dimensional topological 

insulators 
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To guide the search for the Majorana fermion, we theoretically study superconduc- 
tor/topological/superconductor (S/TI/S) junctions in an experimentally relevant regime. We find 
that the striking features present in these systems, including the doubled periodicity of the Andreev 
bound states (ABSs) due to tunneling via Majorana states, can still be present at high electron 
densities. We show that via the inclusion of magnetic layers, this 4-7T periodic ABS can still be 
observed in three-dimensional topological insulators, where finite angle incidence usually results in 
the opening of a gap at zero energy and hence results in a 2ir periodic ABS. Furthermore, we study 
the Josephson junction characteristics and find that the gap size can be controlled and decreased 
by tuning the magnetization direction and amplitude. These findings pave the way for designing 
experiments on S/3DTI/S junctions. 

PACS numbers: 74.78.-w, 74.78. Fk, 74.45. +c, 03.65.Vf 



The prediction of Fu and Kanei that Majorana 
fermions can be realized in superconductor— 3D topo- 
logical insulator structures, boosted theoretical predic- 
tions for the peculiar Majorana fermion properties^—. 
Strong progress has been made in the fabrication of 
two-dimensional 7 - - — and three-dimensional topological 
insulators 7 -^ - — . Recently, the three-dimensional topo- 
logical insulators (TIs) based on the Bi-compounds (e.g. 
Bi 2 Te3, Bi 2 Se3) have already led to the realization of su- 
perconductor(S)/TI/S junctions!*^ and SQUIDs^e. 
From an experimental point of view it is difficult to re- 
alize topological insulator materials with the chemical 
potential at or close to the Dirac point. It is therefore 
highly desirable to have a guiding theory in an experi- 
mentally relevant regime that can pave the way towards 
the verification of the Majorana fermion in S/TI hybrids. 

Here we theoretically study superconductor/three- 
dimcnsional topological insulator Josephson junctions. 
In the calculations of Ref. 1 3.234H it is assumed that 
the Fermi level is close to the Dirac point. In addition, it 
is always assumed that the ferromagnet placed on top of 
a TI, has a magnetization \M\ > /i in these calculations. 
In these systems, an Andreev bound state (ABS) with 
a doubled periodicity is predicted. We consider the ex- 
perimentally relevant regime of high electron densities, 
and show that, despite the chemical potential [i being 
situated far away from the Dirac point, this character- 
istic feature is still present. We furthermore consider 
the presence of a ferromagnetic layer with magnetization 
\M\ < /i on top of the junction, and show that it can 
drastically alter the Josephson characteristics, even when 
(\M\, A) <C /i. This is particularly interesting, since the 
magnetization opens a gap not at the Fermi energy, as 
the superconducting correlations do, but at the Dirac 
point, far away from Ep . We show that a gap in the 
superconducting bound states spectra always opens at a 
finite angle of incidence. However, the size of this gap can 
be tuned and decreased, and can in principle vanish, by 
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FIG. 1: (a) Schematic drawing of a topological insulator (TI) 
with a superconductor (S) and ferromagnet (F) on top. We 
consider the bound states at the surface of the TI with the 
proximity effect from both the superconductor and ferromag- 
net. (b) The energy dispersion at the TI surface (left) and 
at the S side (right). DP indicates the Dirac point. Here the 
case is shown without a F and with an incoming electron at 
the right interface. 



increasing the perpendicular magnetization amplitude. 

The discussion of the bound states of a S/TI/S junc- 
tion in this article is organized as follows: first, we study 
the case without a ferromagnetic layer on top of the TI. 
Then we discuss the bound states with a ferromagnet. 
We will see that a Air periodic Andreev bound state is 
still present in the 3D case but only for one channel. 
This 4-7T periodic ABS is a feature of the presence of Ma- 
jorana fermions. We show the supercurrent obtained by 
the bound states and discuss the observation of a 4-7T pe- 
riodic ABS in a 3D topological insulator. 




FIG. 2: (a) Andreev bound states for different trajectories in a S/TI/S junction. A gap opens at finite angle. The arrows 
indicate which bound state branches are connected to form a 4-zr periodicity. In these branches a Majorana fermion is present 
at (/> = 7r. The legend applies to all the figures. a*t//A = 100, /is/A = 1000 and L/Lo = 0.01 where La = Vfh/firi- (b) 
Air// A = 100,a»s/A = 120 and L/L = 0.01. (c) a*t//A = lOO./xs/A = 120 and L/L = 0.1 and (d) a*t//A = 100,A«s/A = 
120, m z /A = 60.0 and L/Lo = 0.01. (e) Bound state energy for different angles and phase difference <f> = n. Furthermore 
A*s/A = 120, a*t/ /A = 100, m z /A = 0. The energy is oscillating with length. E gap is defined as the distance from E/A = till 
the minimum of the ABS. (f) Bound state energy for fixed phase difference <j> = n and angle using the formula from Ref. l4ll . 



I. THE S/3DTI/S JUNCTION 



the Hamiltonian with a superconducting and magnetic 
proximity effect 



The configuration of the junction we consider, is as 
shown in Fig. [TJa). In the Nambu basis 



H 



H (k) + M A(k) 
-A*(-k) -flS(-k) - M* 



(2) 



where 



H (k) = v F (a x k x + a y ky) - [ij, M — m • er, (3) 

/ at with M the magnetization due to the ferromagnet and 

*= [ip]^h^ v ip]) (1) m = (m x ,m y ,m z ) the exchange field, ai are the Pauli 
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matrices and the index j of the chemical potential is S 
for in the superconducting part and TI in the topologi- 
cal insulator. From this Hamiltonian we calculated the 
eigenvectors of the TI in the presence of the proximimty 
effects. We assume the superconductor to be an s-wave 
superconductor, A(k) = Ae l ^ s , where <p s is the super- 
conducting phase and A = outside the superconductor. 
This Heaviside function of the order parameter simplifies 
the calculations. However, it is not expected that there 
will be a qualitatively difference with the results when a 
self-consistent order parameter is used. As for example 
is shown in Ref. for a d-wave superconductor where 
the induced order parameter is taken into account in the 
normal metal, the Andreev bound state spectrum does 
not change near zero-energy. It is this regime in the spec- 
trum that is relevant for the appearance of a Majorana 
mode. 

The non-superconducting part can be under the influ- 
ence of a ferromagnetic proximity effect. In Majorana de- 
vices the magnetization is taken perpendicular to the TI 
surface, M — m z a z (the resulting eigenvectors are listed 
in the appendix). Magnetization parallel to the interface 
causes only a shift in the wave vector but does not open 
a gap. In the known topological insulators, the Dirac 
point is in the middle of the band gap or close to the va- 
lence band. The chemical potential is usually close to the 
conduction band as the topological insulators based on 
Bi-compounds are not really good insulators yet. There- 
fore, the chemical potential, /iT/,s, is much larger than 
the superconducting gap in experiments. In this case we 
have only normal Andreev reflection at the interface and 
no specular Andreev reflection 3 -^. 

For most proposals it is desired to have a ferromagnet 
on top of the TI (preferably an insulator so that prac- 
tically no current flows through the ferromagnet 
We estimate here how much the gap at the Dirac point 
can be opened by such a ferromagnet. For a magnetic 
moment of nfiB per unit cell of size a 3 , where n is an 
integer and \ib is the Bohr magneton, we can estimate 
the value of the m z a z part of the Hamiltonian. Making 
the assumption that the atoms can be approximated by 
spheres with n elementary dipoles and perfect coupling 
to the TI, we estimate that the opened gap will be about 
0.0002n/a 3 eV where a is in Angstrom. This value is 
typically smaller than the value of the Fermi energy in- 
side the gap (> 0.05 eV). We therefore study the relevant 
regime of m z < hti,s- 

With these assumptions we solve the Andreev and nor- 
mal reflection coefficients at both left and right interfaces 
for incoming electrons and holes by matching the wave 
functions at the interface. Following Kulib 3 ^ (see exam- 
ple in Ref. 33), the wave function in the topological in- 
sulator can be written as 

^ = a^t + + of)- + d^l , (4) 

where the index e and h indicate an electron and hole 
wave function respectively. The coefficients a, 6, c, d are 



related by 

ty+e*^ 1 = r eh , r ot + e t|ke|£ +r Wl!r # h -e-l k ^, 
^+ e -iM£ = reejcij-e^+rhejbip+e-^ 1 , 

where L = ^ cos 6. The second indices of the wave func- 
tions refer to the right (r) and left (1) interface. The right 
interface is placed at L/2 and the left interface at —L/2. 
9 is the angle of incidence from the TI to the S where 
zero angle means orthogonal to the interface. The wave 
vectors kh and k e are the wave vectors of the hole and 
electron respectively. They are given by 

\k h \ = (fj, TI - Ef - ml/vp, 

|k e | = ^(uti + E) 2 -ml/v F , (6) 

where vp is the Fermi velocity. The mismatch between 
these wave vectors and the wave vector in the supercon- 
ductor causes an effective barrier at the interface. The 
superconducting gap can be neglected in this mismatch 
as A < fJ-Ti.s- Conservation of (due to translational 
invariance) then gives for the angle of transmission in 

the superconductor: 9s = arcsin ^s'm9y // ^ I — ml/ ' ns^j ■ 

Solving Eq. (0) gives the energy as a function of the 
phase difference, </>, between the superconductors. 

II. ANDREEV BOUND STATES 

First, we consider the case with no ferromagnet on top 
of the TI. For perpendicular incidence we find a Air peri- 
odic ABS with a gapless dispersion, even in the presence 
of a momentum mismatch (the solid curves in Fig. [5] the 
arrows in Fig. [2ja) follow one ABS that is 4-7T periodic). 
However, in the presence of a momentum mismatch, a 
nonzero angle of incidence results in a nonzero scattering 
amplitude and a gap is always present (Fig. [2](a-c)). The 
larger the mismatch between the wave vectors, the larger 
is the gap that opens. For fipi — ^s t ne interface is ef- 
fectively fully transparent and all trajectories give a Air 
periodic ABS. This is a consequence of the model where 
the superconducting gap is neglected. The opening of 
the gap at finite angles is due to finite back scattering 
at non-zero angle of incidence. For a larger mismatch 
between the chemical potentials, the difference in angles 
of the particle in the TI and superconductor is larger. 
This causes also a larger mismatch in the spin direction 
which increases the barrier and hence results in more re- 
flected electrons. Only at zero angle of incidence, back 
scattering is prohibited by the topological nature. So in 
experiments no Att periodicity of the ABSs can be ob- 
tained for all angles; it is a single channel effect as is also 



4 



(a) (b) (c) 




k/n tyn tyn 



FIG. 3: Bound state energy versus the wave vector parallel to the interface for different phase differences. For all the cases a 
nonhelical Majorana quantum wire is seen at cf> = 71". The legend is shown in (a) and holds also for the other two. (a)/ir//A = 
100,ms/A = 1000 and L/L = 0.01. (b) fi T i/A = 100, ^ S /A = 110 and L/L = 0.01 and (c)/x<ri/A = 100, ^ S /A = 
1000, m z /A = 60 and L/L = 0.01. 



concluded by Fu and Kanei for a system with a small 
Hti,s- Note, that even when this 4tt periodic ABS is 
present, it will only be noticeable in AC measurements 
since interactions with the environment already cause the 
system to reside in the lower (2ir) ABS branche s 29 ' 35 " — . 

For a different length (Fig. H](c)) the curve correspond- 
ing to 9 = 2tt/5 is now lower in energy than the bound 
states of 8 — n/4 and 9 = n/3 compared to the graphs 
of Fig. [21a) and (b). In Fig. [5] (e) the bound state en- 
ergy for specific angles at (j) — tt is plotted as function 
of length. We see an oscillating behavior as function of 
length due to a Fabry-Perot resonance. The oscillation 
period is determined by the Fabry-Perot resonant con- 
dition: 2Lkxi cos 6 = 2nn with kxi the wave vector in 
the topological insulator and n an integer—. There is a 
strong similarity to a normal SNS junction. In Fig. [2] 
(f) a plot is made of the bound state energy of a SNS 
junction for a fixed angle of it /A and phase difference of 
7r by determining the pole in the spectral supercurrent in 
Eq. (5) of Ref. EH i.e. 

r„ =0= (i^ fi 2+^)cosh(^) + 

2Kui n n n sinh [ — - — ] — 

(K 2 - 1) nlcos(2k N L) + A 2 cos</>, (7) 

where Q„ = ^/u; 2 + A 2 , u) n = ir (2n + 1) //3, (3 = l/k B T, 
k N = ( M -C/)-k 2 , k s = ^/i -kf,, 

K = fef - v N(s) = *^£21. T„ = corresponds 
with a pole in imaginary space which is equal to the 
energy of the Andreev bound state. The decrease of the 
amplitude is determined by the ratio of the length of the 
junction and the coherence length of the superconductor 
in the topological insulator. It should however be noted 
that varying the junction length will not result in a 
closing of the gap at a certain length. The calculation 
above is valid for every particular angle. When all angles 



are included the oscillations will be averaged out. 

When a ferromagnet is included, the magnetization 
is found to decrease the gap (see Fig. [5] (d)). This can 
be understood by considering the extreme case where 
the magnetization m is close to the chemical potential 
so that the wave vectors of the electrons and holes are 
nearly zero (Eq. (6)). In that case it also follows from 
the conservation of fcii that S is practically zero. Then 
by using the eigenvectors (Eq. (10), (11) and (17)) in 
the appendix and substitute these values of the wave 
vectors and angle into them, the resulting equations 
at the interfaces simplify. From these equations it 
can be seen that there is perfect Andreev reflection. 
Quantitatively it can be understood by noticing that the 
mismatch between the spins for the different particles in 
the system also causes a barrier. By the magnetization 
this mismatch becomes smaller due to the alignment. 



In Fig. [3] (a) and (b) we show the bound states as a 
function of k y , the wave vector parallel to the interface. 
For a phase difference of 7r between the superconduc- 
tors we see that the zero energy mode has a dispersion 
as function of k y . This is also called a nonchiral Ma- 
jorana stated. For large mismatch between ks and hri 
the model of Fu and Kane applies^ where the chemical 
potential is smaller than the superconducting gap. This 
situation resembles the case of a large mismatch at the 
interface in our case as the waves in the superconduc- 
tor are then fully evanescent. For smaller mismatches 
the gap is smaller, which results in smaller slopes in the 
k y — E graphs. A smaller slope indicates a smaller ve- 
locity of this propaga ting mode along the interface as 
already noted in Ref. |42|. From the same reasoning as 
in the previous paragraph, the k y — E graphs have a 
smaller slope if magnetization is included. The result for 
m z /A — 60.0 is shown in Fig. [3] (c). 
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FIG. 4: (a) Normalized Josephson supercurrent as a function of the phase difference between the superconductors. The numbers 
indicated in the inset correspond with the following parameters: (1) fj. 3 /A = 1000,/iti/A = 100,m z /A = 0, L/Lo = 0.01 
,(2)/i s /A = 120,^tz/A = 100,m z /A = 0,L/L o = 0.01 and (3) /i s /A = 120,/i T //A = 100,m z /A = 60.0, L/L = 0.01. (b) 
Dependence of the normalized J c -Rjve/A product on both the magnetization (dashed line) and the chemical potential of the 
superconductor (solid line) separately. [iti/A is kept constant to 100. In the dependence of the magnetization, fis/A is kept 
constant to 120. The temperature is T/T c — 0.01 in (a) and (b). (c) Sketch of an Andreev bound state at non-zero angle 
of incidence. Due to magnetization the gap of the ABS has become so small that through Zener tunneling the electron in 
the lower branch can be promoted to the upper branch around <f> — n ± n2n where n is an integer, (d) The influence of the 
magnetization and fls/A on the value of the gap. For both situations fin / A = 100, 9 = 7r/3 and L/Lq = 0.01. 



III. SUPERCURRENT 

In this section, we numerically calculate the angle- 
averaged supercurrent of the Andreev bound states. We 
consider here only the supercurrents for small junction 
length as longer lengths give no additional features re- 
garding the 47T periodic ABS and the influence of the 
chemical potential and magnetization on this. In the 
work of Ref. HH discretized bound states were used but 
the continuum was missing in the calculation for larger 
length scales. Because we only consider here small length 
scales, only the discretized spectrum has to be considered 

' // » = £ / / ? cosflta " h (if) f T' <8) 

where Jo = eNA/H. Three plots of this normalized 
Josephson supercurrent are shown in Fig. |4] (a). The 
temperature is T/T c = 0.01 with T c the critical temper- 
ature. 



Although a 4n periodicity is present in the Andreev 
bound state for zero angle of incidence, the other chan- 
nels are 2ir periodic. Hence, the 2tt periodic character is 
dominating the angle-averaged supercurrent and there- 
fore this current will be 2ir periodic in measurements. 
Moreover, the thermal equilibrium of the system even 
makes the ABS for zero angle of incidence 2ir periodic 
since inelastic scattering can relax quasiparticles to an 
ABS that is lower in energy^— . This thermal equilib- 
rium is due to the exchange between the bulk supercon- 
ducting electrodes and the Andreev bound state levels in 
the junction^. 

We see that for larger mismatch the supercurrent as 
function of phase has a more sinusoidal shape. For small 
mismatches there is a sharp transition at <f> = tt from 
positive to negative supercurrent which is also the case 
in a normal superconducting junction^. Also for a larger 
mismatch or for a magnetization the slope of the energy- 
phase curves become less steep, causing a smaller super- 
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current for both cases. In Fig. [4] (b) we plot the de- 
pendence of the normalized I c Rm product as a function 
of mismatch in Fermi level and as function of magneti- 
zation. i?Ar is here averaged over all angles. For larger 
mismatch in Fermi level the value of the I c Rng/ A is sat- 
urating towards 0.57T as expected in the tunneling limit. 
For small mismatch the value is ir as is also the case in 
the ballistic limit in normal SNS junctions. The criti- 
cal current is also decreasing for a larger magnetization. 
However, we see that for values larger than m z /A = 60 
the critical current is increasing again. When analyz- 
ing the Andreev bound states we notice that there is 
a competition between the flatting of the bound states 
and the lowering of the barrier both due to magnetiza- 
tion. The latter depends on the relative magnitude of 
the magnetization to the mismatch. In Fig. 2] (b) the 
mismatch of the wave vectors due to a difference in the 
chemical potentials, is relatively small: (Iti/A = 100 
and us/ A = 120. The aligning of the spins for larger 
magnetization can therefore make the interface almost 
transparent. The corresponding Andreev bound states 
also resemble therefore an almost transparent interface: 
small gap and at <f> = and 2ir the energy is E/A = ±1. 
For a larger difference between (is an d (Xti, the effect 
is less and the critical current is monotonically decreas- 
ing for larger magnetization. However, if we analyze the 
normal resistance, the resistance increases for larger mag- 
netization. This is because the spins in the topological 
insulator with the ferromagnet on top are now more mis- 
aligned compared to the topological insulator side with- 
out a ferromagnet on top. The combination of both an 
increasing I c and i?jv, results in a normalized I c Rn&/ 'A 
product of 2ir. 



IV. SIGNATURES OF 4tt PERIODIC ABS BY 
MEANS OF ZENER TUNNELING 

We have seen that the 4w periodic ABS is a single 
channel effect for perpendicular trajectories only. Next 
to it, measuring in thermal equilibrium makes even the 
single Air ABS 2ir periodic because the electrons will fol- 
low the lowest ABS branches, i.e. below E/A = 0. The 
latter can be solved by doing AC measurements such as 
Shapiro step and/or noise measurements. The Aw peri- 
odic ABS will only contribute to the Shapiro steps at a 
voltage equal to nhuj/e, where n is an integer and u) the 
frequency of the applied microwav o 29 i 35 . A 2ir periodic 
ABS will result in Shapiro steps at V = nhui/2e which is 
half the step size of the Air periodic bound state. 

However, due to the presence of just one single Air pe- 
riodic ABS in 3D TIs/superconductor Josephson junc- 
tions out of many, one would expect that a Air periodic 
signature in AC measurements is not visible due to angle- 
averaging. Usually, one can enhance the zero angle con- 
tribution by introducing a physical barrier of finite width 
due to the exponential dependence of the wave function 
on the width. It is however not possible to cancel the 



nonzero angles by introducing a physical barrier between 
the superconductor and TI because of Klein tunneling 
which renders barriers effectively transparent (see for ex- 
ample the discussion of Klein tunneling in graphene in 
Ref. I3TU32I ). So in order to see a Air periodic ABS, the 
ABSs of all angles should be Air periodic. By means of 
Zener tunneling this can be achieved. 

In order to obtain Zener tunneling, a bias voltage 
across the junction is required (which is for example 
the case in Shapiro steps measurements). Due to this 
bias voltage the quasi-particles in the junction can gain 
enough kinetic energy so that they can transfer from a 
lower Andreev bound state to an upper bound state de- 
spite the separation by a gap 45 ' 46 . It is noted in Ref. 
that for large transparency of the interface this can result 
in a Air periodic ABS. So when the gap is small, it can- 
not be distinguished anymore from an Andreev bound 
state without a gap. Hence, when electrons in the lower 
branch of these Andreev bound states with small gaps 
are promoted to the upper branch, it can result in Air pe- 
riodic signatures in AC measurements 4 ^ 4 ^. In Fig. \M, C ) 
a sketch is shown of an ABS with a finite but small gap. 
Due to magnetization this can result in a Air periodic 
signature in for example Shapiro step and/or noise mea- 
surements. Physically the small gap is similar to having 
a finite length in the ID Kitaev model where also a gap is 
present due to the interaction of the Majorana fermions 
at the end a 36 i 50 ~ — . A way to reduce the gaps of all non- 
zero angle of incidence channels in a 3D TI in order to 
enhance the chance of Zener tunneling to get the Air pe- 
riodicity of all ABSs, is by exploiting magnetization as is 
shown in Fig. [IJd). 

The influence of the magnitude of the magnetization 
depends first of all on the relative magnitude of the chem- 
ical potentials to each other as we can see from Fig. |U(d). 
The larger the mismatch the less the influence is of the 
magnetization. Secondly, the influence of the magnetiza- 
tion depends on its magnitude compared to the absolute 
magnitude of the chemical potential. With a larger chem- 
ical potential, the influence of the magnetization is less. 
To get a clearer picture of the influence of the magneti- 
zation, we have plotted the gap in the Andreev bound 
state (at 4> — ir) for several conditions in Fig. [4] (d). 
For an increase of mismatch wc kept the magnetization 
constant to m z /A = 0, length L/Lq = 0.01 and angle 
= 7r/3. A similar result is obtained for other angles. 
Further fin/A = 100 for both graphs. We see that the 
gap is a strong function of magnetization in the beginning 
but saturates at larger mismatches in the chemical po- 
tentials. The graph that shows the influence of the mag- 
netization has a constant (small) mismatch us/ A = 140 
and uti/A = 100. When the magnetization energy is 
half the value of the Fermi level in the TI, the energy 
gap is already decreased by 50% of its original value at 
zero magnetization. However we have estimated before 
that the magnetization energy is typically on the order 
of 1% of the Fermi energy, i.e. fi 3> m z in experiments. 
So only for small mismatches in the chemical potentials 
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it is possible to reduce the gap so that by means of Zener 
tunneling the 4tt periodicity of the ABSs remains for all 
angles. This reduction of the gap can be made visible 
in SQUID experiments as proposed in Ref. [53U&3 . where 
it is shown that the reduction of the gap size results in 
a different critical current modulation of the SQUID as 
function of the applied flux through the ring. These re- 
sults hold even in equilibrium experiments, so that these 
SQUID devices can be used for ABS spectroscopy. 



CONCLUSIONS 



We studied the Andreev bound states and the result- 
ing supercurrent for S/TI/S junctions with and without 
a ferromagnet on top of the TI. In experiments it is of- 
ten the case that /i » A and m z < /i and therefore 
we extended the model (by Fu and Kanei and Linder et 
alM-) towards this regime. The bound states are solved 
by means of the total wave function in the topological 
insulator and its relation to the reflection coefficients, 
providing insight in the process. The important conclu- 
sion is that the results from Fu and Kane^ (e.g. the 
4-7T periodic ABS existing only for 9 = 0) still remain, 
even for large chemical potentials. Therefore, these fea- 
tures are also valid in actual experiments on TIs. The 
Air periodicity of the bound states only remains in a 3D 
topological insulator for zero angle of incidence. This 4n 
feature cannot be observed as all the other angles, which 
give 2ir periodic bound states, causes the 2ir periodicity 
to dominate. However, Zener tunneling can cause a tran- 
sition from the lower branches to the upper branches of 
the Andreev bound states, even when a gap is present. 
We can enhance this Zener tunneling, and hence enhanc- 
ing the 4-7T periodicity of the non-zero angle ABSs, by 
depositing a ferromagnet on top. The magnetization ef- 
fectively lowers the barrier which causes the gap in the 
Andreev bound states to become smaller. 
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Appendix: Eigenvectors and boundary conditions 

The eigenvectors in the topological insulator are cal- 
culated to be 



^2 = P2 



V>3 = P3 



V>4 = n<i 



( -m + y/m 2 + t; 2 |ki| 2 ^ 






\ 



y/m 2 + v 2 \k\ 2 



\ -v|k 2 |e-» J 

' m + \J m 2 + w 2 |k3| 2 \ 
v\k 3 \e l9 


V o / 



—to + y to 2 + t; 2 |k4 1 
v\k 4 \e- t9 



\ 
) 



(9) 



(10) 



(11) 



(12) 



where nj = 1/ ^J2{m 2 + ?j 2 |kj| 2 — m^/m 2 + u 2 |kj| 2 ) 
with j = 1,4 for ipi and ip4 respectively. Further pj = 

l/y / 2(m 2 + w 2 m| 2 + my/m 2 + v 2 |kj| 2 ) with j = 2,3 for 



ip2 and ip3 respectively. The eigenvalues are 



Ei = -n - Vr 
E% = fi — \/m< 



V|ki| 



E 3 



« 2 |k 4 | 



(13) 
(14) 
(15) 
(16) 



■01 and ip3 are the electrons belonging to the lower and 
upper half of the Dirac cone respectively. "02 & n d ip4 are 
the holes corresponding respectively to the lower and up- 
per half of cone. The wavefunction in the superconductor 
is given by 



/ e**y/E-n + v\k.\ \ 



e^e^y/E 



Ae l9 



v\k s 



VE 



fi + v\k s 
A 



(17) 



V ^E-n + v\k s \ j 



with an energy given by E = -\/A 2 + (w|k s | — /i) 2 . 

The reflection and transmission coefficients can differ 
at both interfaces due to magnetization. If we consider 
the electrons at the upper cone as depicted in Fig. Q] 
(a) we need the wave functions ip2 and 1P3. Taking the 
direction of the particles into account in the angle 9 in 
the TI and 9 S in the superconductor we arrive at the 
following set of equations: 
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Right interface, incoming electron 



Pm,3 (1 + Tee) 
( id -%B\ 

Pk,3 (e - r ee e ) 

Pm^eh 



-i9 

~Pk,2e r eh 



oi<l> 



Left interface, incoming electron 



(Citee + diteh) ■ 



\cit ee e l 3 - dit eh e~ 



2VE 

e i<t> 

A 



2-s/E V c 2 
A ( 1_ 1_ 

2-s/E \ c 2 d 2 J 



e 3 e 

tee ~t~ t e h 



-i0„ 



(18) 



Pm,3 (1 + r ee ) 
Pfc,3 (~e- 4e + r ee e 10 ) 



(dlt e h + Cit ee ) . 



(dit eh e ie ° - Cl t ee e- ie °) 



where 



Pm,2 r eh 



|k se | - n/v + Ve 2 ~ a 2 /v, 

\k sh \ = p/v-y/E*-*?/v, 



Pk,2 r eh.e 



it) 



2\/E 
2VE 

A / e l ° 3 e-' w 
2\fE V d 2 c 2 

A / 1 1 

teh~r + 



i "en , i "ee 

V d 2 c 2 



(21) 



Cl 


= 




- v|k se |, 




di 


= y/E 




- v|k sh |, 




C2 


= Ve 




- u|k se |, 


Left interface, incoming hole 



d 2 

Pm,j 

Pkj 



y/E - n + v\k sh \, 



m + y/m? + v 2 |kj| 2 



^2(m 2 + v 2 |kj | 2 + ra^Jm 2 + v 2 |k,-| 2 ) 



2(m 2 + v 2 \kj | 2 + my/m 2 + i> 2 |kj| 2 ) 



.(19) 



Right interface and incoming hole 



Pm,3The 

-Pk,nr he e~ te 

Pm,2 (1 + Thh) 
( i6 -i6\ 

Pk,2 (e - r hh e j 



(Cit ee + dlteh) ■ 



( Cl t ee e ie ° - d^ehe- 103 ) , 



e i<j> 
2\/E 

e i<f> 

2\/E 

A / e i9s e~ id - 

7= 1 t e e ~\~ teh ; 

2\fE V c 2 d 2 
A / 1_ l_ 

2y/E V C 2 ^2 



Pk,3rh e e 

Pm,2 (1 + rfc/i) 

Pfe,2 (-e" ifl + r hh e ie ) 



r— (dlteh + Citee) i 

2v £/ 

(dit eh e <9 ' - ateee- ie °) , 

2v ii 

A / e ie * e~ l9 
2VS V <fc c 2 



-4= (teh^-+t ee - )• !22 i 

2V-B V "2 ci 



(20) 



These equations can be used to calculate the coefficients 
that are used in Eq. ([5]) in the main part of the paper. 
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